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OBJECTIVE 


Determine  the  feasibility  of  sealing  photopie  almospherie  transmission  to  the  infrared 
wavelengths  for  a marine  environment.  Speeifieally.  eonduet  measurements  of  marine 
aerosol-si/e  distributions  and.  by  using  Mie  theory,  derive  sealing  laws  to  relate  the  photopie 
and  infrared  e.xtinefion  eoeffieients  for  differing  visibility  eonditions. 


RESULTS 


1 .  Simple  linear  relationships  ean  be  derived,  by  using  Mie  theory  and  measured 
aerosol-size  distributions,  which  relate  the  photopie  and  infrared  atmospherie  extinetion 
eoeffieients. 


2.  Similarly,  linear  relationships  ean  be  derived  whieh  relate  the  wavelength-dependent 
extinetions  to  the  atmospherie  liquid  water  eontent. 

3.  Different  linear  relationships  foreaeh  wavelength  are  derived  depending  upon 
whether  the  meteorologieal  visibility  is  less  than  or  greater  than  about  4 kilometres. 

4.  Measured  meteorologieal  visibilities  are  elo.sely  related  to  measured  almospherie 
liquid  water  eontents.  1 his  suggests  that  simple  visiometer  or  laser  baek.seatter  ( lidar)  measure- 
ments ean  be  used  with  sealing  laws  to  estimate  atmospherie  transmittanee  for  infrared 
wavelengths. 


Compari.son  of  extinetions  predieted  analytically  on  the  basis  of  accepted  fog 
aerosol  models  and  those  numerically  derived  from  measured  distributions  indicates  that  the 
eontributions  from  particles  having  diameters  greater  than  .10 /am  may  signifieantiv  alter  the 
derived  sealing  laws. 


RECOMMENDATIONS 


I.  ( ontiiuie  measurements  of  aerosol-size  distributions  in  a variety  of  marine  environ- 
ments. l-.xtend  the  measurement  capabilities  of  the  Knollenberg  spectrometer  to  include  those 
particles  having  diameters  exceeding  .10  urn. 


-.  Obtain  airborne  measurements  of  aerosol-size  distributions  sueh  that  sealing  laws 
ean  be  obtained  whieh  include  horizontal  and  vertical  variations  in  aerc^sol  concentrations. 

.1.  Extend  the  instrumentation  capabilities  to  include  both  tnilsc-d  and  EM->  w li.l  .rs  to 
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INTRODUCTION 


The  perfoniiance  of  existing  and  planned  Navy  shipboard  electro-optical  (EOl 
systems  which  utilize  infrared  (IR)  wavelengths  is  severely  limited  by  the  constituents  of  the 
atmosphere  and  the  inhomogeneitics  in  their  concentrations.  The  ability  of  an  EO  system  to 
perform  a function  ie.  detect  or  identify  a target  within  specified  probability  limits  depends 
upon  the  target  irradiance  signal- to-noise  (background  plus  internal)  ratio.  The  irradiance 
signal-to-noise  ratio  is  degraded  by  atmospheric  extinction,  which  is  the  sum  of  absorption 
and  scattering  by  aerosols  and  molecular  gases.  The  contributions  to  the  extinction  from 
molecular  water  vapor  are  important  at  IR  wavelengths  and  can  usually  be  assessed  by  stan- 
dard meteorological  measurements.  More  difficult  to  assess  are  the  contributions  from  aerosols 
and  changes  in  their  size  distributions,  which  depend  to  a large  extent  upon  relative  humidity 
changes  In  situ  measurements  of  aerosol-size  distributions  made  with  groundbased  or  air- 
borne devices  ( Knollenberg  spectrometers,  impactors,  etc),  together  with  Mie  theory,  provide 
a way  for  estimating  atmospheric  aerosol  extinction  at  ditterent  wavelengths.  These  devices 
are  usually  expensive  and  impractical  for  operational  licet  uses.  Needed  to  assist  system 
planners  and  users  in  predicting  the  performance  of  shipboard  EO  systems  are  ( 1 ) remote 
sensing  techniiiues  for  relating  meteorological  parameters  (which  are  unaffected  by  the  local 
shipboard  environment)  to  the  extinction  at  visible  wavelengths,  and  (2)  scaling  laws  which 
then  relate  the  visible  extinction  to  that  in  the  IR  band. 

Bieberman  ct  al  ( l‘)77)  recently  published  results  based  on  continental  aerosol-siz.e 
distribution  measurements  at  Crafenwohr.  l ederal  Republic  of  Cermany.  during  the  1975-76 
winter,  relating  the  liquid  water  content  to  the  aerosol  extinction  coefficient.  They  suggest 
that  a scaling  of  the  photopic  transmission  to  the  IR  wavelengths  is  possible  and  that  the 
scaling  laws  are  independent  of  the  particle-size  distributions.  Erom  the  earlier  work  of 
Houghton  ( 1939)  and  others  as  referenced  by  1 Idridge  ( I9(,(,),  it  is  recognized  that  a gen- 
eral inverse  relationship  exists  between  the  atmospheric  liquid  water  content  and  the 
visibility  as  defined  by  Koschmieder's  formula  ( = 3.92  where  is  the  aerosol 
extinction  coefficient).  However,  the  relation  between  the  visibility  or  extinction  coetti- 
cient  and  the  liquid  water  content  is  inherently  dependent  upon  the  aerosol-size  spectrum  tor 
clear  weather  conditions.  Only  for  the  poorer  visibility  conditions,  such  as  fog,  is  the  relation- 
ship determined  primarily  by  the  largest  droplets  present  and  not  so  much  the  individual 
distribution. 

In  this  report,  aerosol-siz.e  distributions  measured  prior  to  and  during  a coastal  marine 
fog  are  presented.  Together  with  Mie  theor>’.  the  measured  distributions  are  used  to  derive 
scaling  laws  which  a-late  the  liquid  water  content  to  (he  photopic  and  IR  extinction  coeffi- 
cients for  the  differing  visibility  conditions. 


MEASUREMENTS 

In  addition  to  the  standard  surface  meteorological  observations,  aerosol-siz.e  distribu- 
tions are  routinely  made  at  the  NOSC'  remote  sensing  facility  located  on  the  sea  side  of  Point 
Loma  in  San  Diego.  California.  The  frequently  occurring  coastal  marine  fogs  related  to  Santa 
Ana  or  stratus  lowering  events  provide  the  opportunity  to  make  atmospheric  measurements 


Biberman.  LM.  Roberts,  RK.  and  Scekamp,  LN,  A comparison  of  electro-optical  technologies  for  target 
acquisition  and  guidance.  Institute  for  Defense  Analysis,  paper  P-1218.  January  1977 
Houghton,  HG,  On  the  relationship  between  visibility  and  the  constitution  of  clouds  and  fog.  Jour  Aero 
Sci,6,408,  1939 

Eldridge,  RG.  Haze  and  fog  aerosol  distributions.  Jour  Atmos  Sci,  23.  oO.t,  1966 


in  a variety  of  rediiccd  visibility  conditions.  For  tliis  analysis  aerosol-size  distributions  are 
used  which  were  obtained  during  a 10-hour  period  on  I.S  February  F)77  when  the  visibility 
varied  from  approximately  0.1  to  6 km.  The  distributions  were  measured  by  using  an 
aspirated  Knollenberg  ASSP-lOO  spectrometer  probe  mounted  atop  a building  approximately 
45  metres  above  mean  sea  level.  The  spectrometer  si/.es  particles  in  four  overlapping  bands 
which  are  individually  divided  into  15  different  diameter  inter\'als.  Spectra  are  measured  in 
each  range-band  every  second  and  digitally  recorded  on  magnetic  tape.  Fhe  digitized  data 
are  subsequently  processed  on  a Data  General  C’orimration  NOVA  XOO  computer.  Individual 
spectra  in  each  range-band  were  averaged  for  4 minutes,  and  then  combined  to  give  a diam- 
eter coverage  from  0.45  to  2^). 4 /am.  The  liquid  water  density  was  also  computed  by 
integrating  the  averaged  size  distributions  assuming  the  aerosols  to  be  predominately  spherical 
water  droplets.  Throughout  the  10-hour  period  the  meteorological  visibility  was  measured 
with  a forward  scatter  fog  visiometer  (Meteorology  Research  Inc  model  I 5X0)  which  was 
located  adjacent  to  the  Knollenberg  spectrometer.  Figure  I shows  the  liquid  water  density 
computed  from  the  measured  size  distributions  at  5-minute  intervals  for  the  lO-hour  period 
and  the  measured  meteorological  visibility.  Good  inverse  correlation  is  evident  between  the 
liquid  water  content  and  the  visibility  for  the  entire  period. 

F.xamples  of  the  measured  size  distributions  are  shown  in  figure  2.  IVior  to  (02  1 5/. 
and  0245Z)  and  after  (O.^OOZ  and  03  I 5Z)  the  advection  of  the  tog  over  the  recording 
station  the  distributions  approximate  inverse  power  laws  for  diameters  greater  than  1-2  /am 
with  diameter exponentsof  3.b5.  3.55.  2.53.  and  2.73.  respectively,  fhe  number  densities  at 
diameters  greater  than  5 /am  increased  tw'o  orders  of  magnitude  tiuring  the  log. 


DATA  ANALYSIS 

A computer  program  (compiled  by  Stone.  Ib7(i)  was  used  to  calculate  the  extinction, 
scattering,  and  absorption  coefllcients  for  optical  propagation  through  aerosols  and  at  wave- 
lengths of  0.53.  l.()(i.  3.5.  and  I0.(i/am.  The  contributions  to  absorption  by  molecular  gases 
have  not  been  included  in  the  calculations.  Basically,  the  program  evaluates  the  integral 
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(r/X.  rj)  dr 
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where  the  coeflicients  of  concern  are  designated  by  |3  and  the  subscripts  ae  and  as  reler  to 
aerosol  extinction  and  scattering,  respectively.  I'he  aerosol  absorption  coefficient  is  then 
determined  by 


/^aa=^ae-^as 

In  equation  ( I ).  n(r)  is  the  number  of  particles  per  unit  volume  having  radii  between  r and 
r + dr.  K^.(r/X.  rj)  is  the  total  cross  section  for  extinction  or  scattering  determined  by  Mie 
s 

theory  and  normalized  to  the  spherical  particle  geometrical  cross  section,  rj  is  the  complex 
index  of  refraction  of  liquid  water  at  a wavelength  X as  presented  by  Deirmendjian  ( 19(j4). 

Slone,  WR,  MIF  and  SCATCO:  Programs  to  compute  the  scattering,  extinction,  and  absorption  properties 
of  optical  propagation  througli  particle  distributions.  MlitiATliK  Corporation.  San  Diego.  CA.  Report 
R 2005 -070- 1 - 1 . I November  F)76 

Deirmendjian.  D.  Scattering  and  Polarization  Properties  of  Water  Clouds  and  Ha/es  in  the  Visible  and 
Infrared.  Applied  Optics,  .t(2).  187.  t')b4 

-t 


Actually,  equation  ( 1 ) is  integrated  hy  using  Simpson’s  rule  over  une(|ually  spaced  intervals. 

In  the  integration  a third-order  polynomial  lit  is  used  to  interpolate  the  measured  distrihu- 
tions  between  the  diameter  interval  midpoints  to  tit  the  size  parameters  for  which  the  Mie 
coefficients  are  given. 

Figure  3 shows  the  extinction  coefficients  calculated  at  lO-minutc  intervals  for  the 
initial  4-hour  period.  The  extinction  coefficients  closely  follow  the  liquid  water  density 
variations,  as  is  to  be  expected,  since  both  are  determined  by  integrating  the  size  distribu- 
tions. Of  significance  are  the  different  relations  between  the  visible  (0..‘s3  /am)  extinction 
and  the  extinctions  for  the  IR  wavelengths  before  and  after  the  onset  of  fog.  Bcibre  the  fog 
the  O..S3-/am  wavelength  suffers  the  greatest  extinction,  while  after  the  fog  onset  there  is 
little  difference  between  the  extinctions  of  the  three  smaller  wavelengths.  The  ratio  of  the 
0..33-  and  10.6-^im  extinctions  changes  by  approximately  a factor  of  after  the  onset  of  tog. 
riiese  differences  illustrate  the  visibility  dependence  on  any  scaling  laws  one  would  derive 
from  the  measured  size  distributions.  This  is  further  evidenced  in  figure  4,  where  the  ratios 
of  the  scattering  to  absorption  coefficients  at  3.5  and  lO.O  /Jin  are  shown  as  a function  of 
time.  Before  and  after  the  fog  onset  the  extinction  at  3.5  /am  is  primarily  absorption  while 
at  1 0.(1  /am  absorption  dominates  before  the  fog  and  scattering  dominates  after  the  tog  onset. 
However,  these  results  may  be  significantly  altered,  particularly  for  the  t0.(i-/am  wavelength, 
when  contributions  to  the  size  distributions  frvim  droplets  with  diameters  greater  than  30  /am 
are  included  in  the  calculations. 

In  figure  5 the  aerosol  extinction  coefficients  at  3.5  and  10. (>  /am  are  plotted  versus 
those  at  0.53  ao'i.  I he  1 .0(i-/am  extinctions  are  excluded  because  they  ditfer  only  slightly 
from  those  for  0.53  /am.  Regression  analyses  of  the  data  show  different  linear  relationships 
between  the  photopic  and  IR  extinctions  for  visibilities  less  than  and  greater  than  4 km. 
Similarly  in  figure  (i.  different  linear  relationships  are  derived  between  the  liipiid  water 
density  and  the  0.53-/am  extinction  for  water  densities  less  than  and  greater  than  0.005  g m \ 
.Mthough  the  correlation  coefficients  ( r)  for  the  different  regressions  are  good,  there  is  con- 
siderable scatter  in  the  data  before  the  fog,  especially  for  the  I0.(i-pm  wavelength,  which 
again  illustrates  the  dependency  on  the  individual  distributions,  fhis  dependency  can  be 
removed,  somewhat,  if  the  data  in  tigures  5 and  (i  are  normalized  to  the  total  number  of 
particles  per  unit  volume.  N.  in  each  ilistribution.  as  shown  in  figures  7 and  X.  respectively . 
However,  to  he  of  practical  value,  any  scaling  laws  derived  from  the  figures  would  retpiire  as 
inputs  the  individual  size  distributions  as  well  as  the  liquid  water  content. 

From  the  regression  analy  ses  in  figures  5 and  (i,  the  following  eipiations  can  be 
derived  which  relate  the  aerosol  extinction  at  each  wavelength  to  the  liquid  water  contents 
\V|  > 0.005  g nr^  > W V. 
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Figure  3,  Aerosol  exiinclion  coelllcieni  and  Ijrjuid  water 
content  as  a function  of  time. 
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III  figure  'i,  the  extinctions  derived  for  different  water  densities  by  using  the  15  l-'ebruary 
scaling  laws  are  compared  witli  those  calculated  from  selected  individual  si/e  distributions 
measured  during  different  fog  periods  which  occurred  8 and  I 7 l-ebruary  |b77.  The  good 
aga'ement  indicates  a consistency  in  the  size  distributions  for  the  different  fog  periods  and 
demonstrates  the  utility  of  the  scaling  laws.  It  is  evident  that  IR  extinction  coeftlcicnts  can 
lie  estimated  from  simple  scaling  laws  using  visiometer  measurements  of  the  photopic 
extinction  or  remote  sensing  techniques  such  as  laser  ranging  devices  (lidar).  as  suggested  by 
Roberts!  I‘)76).  whose  backscatter  profile  is  sensitive  to  the  liquid  water  conlcnl  along  the 
propagation  path. 

COMPARISON  WITH  THEORY 

I here  is  little  dillerence  between  the  ().5.T  and  .T5-/L:m  extinctions  shown  in  ligure 
for  water  contents  greater  than  about  0.0.^  g m-\  1 he  fallotf  in  the  curves  at  lO.b  /am  mav 
simply  be  caused  by  limiting  the  maximum  radius  to  14.2  |im  ( the  midpoint  of  the  highest 
range-band  of  the  Knollenberg  spectrometer).  Had  larger  drop  sizes  been  measurable  by  the 
Knollenberg  probe,  their  contribution  to  the  IO.()-;am  extinctions,  due  to  increased  absorp- 
tion. might  have  llattened  the  curves  over  the  wavelength  band  considered  for  the  higher 
water  densities.  However,  the  near  flatness  ot  curves  between  0.5.^  and  .T5  pm  allows  an 
analytical  expression  for  extinction  to  be  obtained  which  is  primarilv  ilependent  on  the 
liquid  water  density  and  not  so  much  i>n  'he  individual  size  distributions  for  heavy  fogs. 

On  the  basis  of  experimental  results  obtained  worldwide.  Junge  ( l‘)52)  concluded 
the  size  distributions  of  "large  and  giant"  particles  could  be  approximated  by  the  simple  law 

n(  r)  = ('  r'f  (4) 

where  ('  and  p are  constants.  I iiuation  I I ) can  now  be  written 


K( r X ) dr 


Changing  the  variable  of  integration  to  a = r|  X one  obtains 

a~~f  K(a)  da 


P j;- 

Jo 


where  the  lower  limit  of  integration  is  ignored  because  r^,  r| . Now  |j  = KX,  a)  and 
a = f(r|.  X)  such  that  the  total  differential  of^l  is  defined  by 


)P  <10  da 
OX  ^ Oa  OX 


dX 


(O) 
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L 


^ae  using  scaling  taws  derived 

from  15  Feb  drop-size  distributions 

Calculated  from  drop-size  distributions 
measured  on  8 Feb  ( • ) and  17  Feb  (XI 


Performing  tlie  requircil  operations  on  (h)  ami  substituting  into  ( 7 ),  Kurniek  el  al  ( I ‘>(>01 
arrived  at  the  following  expression  (which  is  similar  to  one  derived  by  Junge): 


d(Cn(?y^^.) 
d(Cn  X) 


i'lr  r 


3-1 


(3-p) 


/^ae 


K(r] 'X) 


(XI 


By  using  (4),  the  constant  C in  (X)  can  be  determined  from  the  licjuid  water  density  given  b> 


W = C 


4 3 TT  r'^  P p dr 


(‘>1 


where  p is  the  water  density.  For  p < 4,  and  ignoring  the  lower  limit,  integration  of  ( ‘>1  yields 

O.V.sttC  r|“P 

W ^ (10) 

(4-p) p 


from  which  C can  be  obtained,  and  (X)  then  becomes 


d(tn  .7.^<4-p)  K(r|  X)  \V 

^ =(3-p)-  . (11) 

d(VnX)  '"iPP’.ie 

Junge  ( I 952)  deduced  that  p 2 for  fogs.  .Mso.  when  r j X ^ 1 . K(  T]  X ) ^2.  F igure  X sliov^s 
that  d(Cn  l3..,.)/d( Kn  X)  ^ 0 for  the  higher  water  densities  such  that  (II)  can  be  solved  lor  the 

clC 

extinction  coefficient  to  give 


==  21  1 W.  kin~'  (12) 

where  was  taken  as  14.2  pm  and  VV  has  units  of  g nr\ 

In  figure  10  the  extinctions  derived  by  using  the  Junge  fog  model  are  compared  will) 
the  regression  analy.ses  (equation  (3))  for  the  higher  water  densities.  For  the  measured  si/e 
distributions  the  condition  that  rj/X^  I is  satisfied  only  for  the  0.33- and  l .0(>-pm  wave- 
lengths. The  close  agreement  between  the  3.5-pm  regressions  and  the  analytical  values  is 
probably  fortuitous  since  K(r|/X)  for  this  wavelength  is  1.34  rather  than  the  asymptotic 
value  of  2.  The  10.6-pm  regressions  (not  shown)  differ  from  the  analytical  values  by  an  order 
of  magnitude.  There  is  fair  agreement  between  the  ()..33-  and  I .()(>-pm  regressions  and  the 
analytical  values  for  10““  < W < 10“' . fhe  differing  slopes  most  probably  result  from 
neglecting  the  contributions  from  larger  particles  to  the  scaling  laws.  Whether  or  not  particles 
with  radii  greater  than  14.2  pm  existed  in  the  fog  analyzed  here  cannot  be  determined.  How- 
ever. Junge  ( I ‘>.“>2)  states  that  r|  tor  Fiuropean  fogs  is  about  20  pm  and  Houghton  ( l‘>3‘>) 
reports  values  exceeding  50  pm  for  coastal  and  mountain  fogs  in  the  eastern  United  Slates 
Clearly,  while  these  data  demonstrate  the  utility  of  the  scaling  laws,  additional  measurements 
need  to  be  made  in  a variety  of  different  marine  environments  to  include  in  the  derived 
sealing  laws  those  particles  with  diameters  exceeding  30  pm.  In  addition,  airborne 

Kurniek,  SW,  Zittcr,  RN,  and  Williams,  l)B.  Allcniiation  of  infrared  radiation  by  Fogs,  Jour  Opt  Soe  .\in. 

50,  578,  1960 


measurements  of  aerosol-size  distributions  need  be  made  such  that  scaling  laws  can  be 
obtained  which  include  horizontal  and  vertical  variations  in  aerosol  concentrations. 


J 


I 


figure  10.  Comparison  of  extinction  coefficients  determined 
frirm  regression  analyses  and  those  calculated  by  using  tltc 
■lunge  fog  drop-si/e  distribution  model. 
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